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Abstract 

We calculate the branching ratio of the non-resonant B —^ n“7r+7r+7r“ 
decay using a simple model based on the framework of the factoriza¬ 
tion approach where one can ommit the nonfactorization effects. In 
naive factorization approach, there are only two tree diagrams for this 
decay mode. In the first diagram, the matrix element of decay mode 
is factorized into a, B D form factor multiplied by a Stt decay con¬ 
stant and in the second diagram, the matrix element is factorized into 
a B —>■ Dtt form factor multiplied by a 27r decay constant. We assume 
that in the rest frame of B meson, the D meson remains stationary, 
so we obtain the value (3.47 ± 0.14) x 10“^ for the branching ratio of 
the B —>■ n“7r+7r+7r“ decay mode, while the experimental results are 
(3.9 ± 1.9) X 10-3. 

1 Introduction 

The three-body and four-body decays of B meson are more complicated than 
two-body case, specifically where both non-resonant and resonant contribu¬ 
tion exist. In this types the theoretical momentums of the output particles 
are not directly calculable. In four-body cases, the momentums and form 
factors are written in terms of five variables: two invariant-mass of particles 
variable and three angles. 

In the previous works [HE], for some specific three-body decays assumed 
that one of the output mesons due to its heavy weigh, remains stationary 
and the two other output mesons move back to back and the amplitude 
of this three decay modes can be calculate. In this work we extended this 
assumption to non-resonant four-body decays and calculate the branching 
ratio of R —7- D-7r'*“7r'*“7r- decay in naive factorization approach where we 
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can ignore the non-factorization effects. We assume that the D output me¬ 
son, due to its heavy weight, remains stationary and the three pion mesons 
move in a plane and we have directly resolved their momentums and we cal¬ 
culated the branching ratio of D~Tr~ and predicted (3.47±0.14) x 10“^ 
which is compatible with the experimental result (3.9 ± 1.9) x 10“^. 

This paper is organized as follows. In section 2 we consider the amplitudes 
of the Zl~7r'*“7r'*“7r“ decay and we present the transition matrix needed and 
the kinematics of the four-body decays. In section 3 we give the numerical 
results, and in the last section we have a summary. 

2 Non-resonant amplitudes of the B —> 

For the four-body B —?■ ZI“7r'*'7r^7r“, the Feynman diagrams are shown 
in Fig. 1, where both panel show the b —)• cud tree transitions Elia [5]. 
Under the factorization approach, the contribution of the diagram (a) in 
decay amplitude consists of the < i? —)• D~ > x < 0 —)• Trvrvr > distinct 
factorization term, where < B ^ D > denotes the form factor of B to 
D meson transitional, and the contribution of the diagram (b) consists of a 
< i? —> TTTT > X < 0 —)• TTTT > distiuct factorization term, where < i? —)• tttt > 
denotes the two-meson transition matrix element. So the matrix elements 
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Figure 1: Feynman diagram for the B ^ D Tr~^7r~^7r decay. 
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of the B ^ D decay amplitudes are given by 

\M\ = i^^VcbVud\ai < T^'T^A{du)-\^ >< D\{ch)-\B > (1) 

+02 < Z?7r|(c5)_|0 >< 'n'K\{du)-\B > |, 

Where in factorization, the general expression of the Oj quantities in terms 
of effective Wilson coefficients is [7] 

1 

Ol = Cl + -C 2 , 

1 

02 = C2 + -Cl. (2) 


2.1 Form factors and decay constants 


According to Fig. 1(a), the form factor of B meson to D meson is defined 
as fallow [6] 


< D{p4)\{cb)^v-A)\B{P) > = {PB +P4)fi - ^ -P4)ii Fi(g^) 


9 2 


{PB - PiY 

and the point like 3-body matrix element is assumed to be [8] 


(3) 


< TT{pi)7r{p2)TT{p3)\{du)(^V-A)\0 > 


2i 


P2^l 


{PB -P4.)-P2 
{PB - PaY - ml 


{PB - P4)^, 


(4) 


According to Fig. 1(b), for the current-induced process, the two-meson 
transition matrix element < TnT\{du)^v-A)\B > hos the fallowing general 
expression BM 


< Tnr\{du)(^V_A)\B > = ir{pB - Pi - P 2 )tM + i^+iP 2 + Pi)iJ. + ii^-iP 2 - Pi)tJ., 

(5) 

where the form factors r,uj^ and iO- for the non-resonant decay are given 
by 

_ Ib Ib PB-{P 2 -Pl) I niB {PB-Pl)-Pl 

2/2 /2 ■ (p^ - Pi - P 2 Y -m% fn mB* ■ {pB -PlY - 

4^9‘^fB mBmBi Pi-P 2 -pi-{pB-pi)P 2 -{pB-pi)/m\,^ 

^ ' {PB - Pi - P2Y - m%' (ps - Fi )2 - m|. 

( 6 ) 
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fn' {PB-Plf - rUB* 


(7) 
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_ g fB*'mB*y/mBmB* 

{PB - Pi)^ - rriB* 

The tow-body < Z) 7 r|(c 6 )(y _^)|0 > matrix element can be related to the D 
and TT matrix element of the weak interaction current: 

< D{pi)TT{p^)\{ch)(y-A)\^> = (pa-P 3)F^''{q‘^), (9) 

where is the D to vr transition form factor and determined from 

experiment m- 


1 + 


{PB-Pl)-Pl 


m 




( 8 ) 


2.2 The kinematics of four-body decay 

For four-body B decays, we can work in the rest frame of the B meson. By 
following the definition given in mm, when B meson decays to hnal-state 
mesons, the kinematics of the four-body decays can be described in term of 
five independent variables where the momenta and masses are defined as 


P{p,M) -A Pi{pi,mi)P2{p2,rn2)P3{P3,rn3)P4{p4,m4), (10) 

with p^ = and pf = rrif where i = 1,2,3 and 4. We choose the set 
of independent variables as {si 2 , 534 , 0 i, ^ 3 , (/>} which have the following ge¬ 
ometrical meaning (see Fig. (2)) 

• S 12 = {pi +P 2)^5 is the invariant-mass of particles 1 and 2 ; 

• S 34 = {p 3 +Pa)‘^, is the invariant-mass of particles 3 and 4; 

• ^ 1 (^ 3 )) is the angle between the three-momentum of particle 1 (particle 
3) with respect to the direction of pi 2 = pi +P 2 {psa = Ps +Pa) defined 
in the rest frame of the decaying particle (see fig.(2)); 

• 4> is the angle between the planes defined by particles (1, 2) and (3, 4) 
also in the rest frame of the decaying particle. 


With this choice of kinematics, the differential decay rate in the rest frame 
of the decaying particle can be written as: 


dT 


^dsi 2 ds 34 d cos Old cos O^dip, 
4 ( 47 r)°m^ n! 


( 11 ) 


where / 3 i 2 (/ 334 ) is the velocity of particle 1 (particle 3) in the center of mass 
frame of particles 1 and 2 (3 and 4) and X = [(P^ —S '12 — 5 ' 34 )^ — 45 i 25 ' 34 ]^/^, 
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Figure 2: Kinematics of four-body decays in the rest frame of the decaying 
particle, = 0. We have defined pij = Pi+Pj, such that pi2 +P34 = 0. 


the n! factor accounts for identical particles in the phase space and \Ai\ is 
the amplitude of the decay (Eq.(l)). 

For considering the B —)• D^vr+Tr+vr” decay mode, we work in the rest frame 
of the B meson and assume that the D~ meson due to its heavy weight, has 
a small momentum and can easily be neglected, so the three pions move in a 
plane. Hence the momentums of B meson and final-state mesons are given 
by 

Pi = (Pnl^Pnl), P4 = {mD,0), 

P 2 = (P°2>P7r2), PB = {ms, 0) 

P3 = (p°3,P7r3), 

where = Pt ,2 = PnS and ph = p °2 = pIs = - mn)/^- Now we can 

easily integrate the decay rate with respect to all independent variables and 
we obtain a simple experiment 


F 


Xj3i2(33A 

4(47r)®m^ 


|A4|2-^S12S34. 


( 12 ) 


3 Numerical results 

The Wilson coefficients Ci have been calculated in different schemes. In 
this paper we use consistently the naive dimensional regularization (NDR) 
scheme. The values of c* at the scales p = at the leading order (LO) are 

m 


Cl = 1.117 C2 = -0.268. 
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The parameter g in the form factors, determined from the D* —)• Dtt decay 
is 


g = 0.3 ±0.1. (13) 

For the elements of the CKM matrix, we use the values of the Wolfenstein 
parameters and obtain |14| 


Kd = 0.9745 ^^6 = 0.0415 

The meson masses, decay constants and form factors needed in our calcula¬ 
tions are taken as (in units of GeV) [T5l 116] 

mB = 5.279 ms* = 5.413 rriB = 1.87 m-,^ = 0.139 

/b = 0.176 /b/ = 0.220 /, = 0.13 = 0.52 

Ff = 0.52. 

Using the parameters relevant for the B —)• F“7r+7r'*'7r“ decay, we calculate 
the branching ratio of this decay which is 

B{B F"7r+7r+7r“) = (3.47 ± 0.14) x 10"^ (14) 

where the experimental value is |15j 

B{B F"7r+7r+7r“)expt = (3.9 ± 1.9) x 10"^. (15) 

4 Summary 

In this work, we have considered four-body B —?■ F“7r'’'7r'*'7r“ decays under 
the assumption that the D meson has small and negligible momentum due 
to its heavy weight. This decay mode have contribution of ai and a 2 coef¬ 
ficients. According to the Feynman diagrams, each contribution factorized 
into a form factor and a decay constant. Under the assumption of negligible 
momentum of D meson, the three output pion mesons move in a plane and 
we have directly resolved their momentums and we calculated the branching 
ratio of and predicted (3.47 ±0.14) x 10“^ which is compatible 

with the experimental result (3.9 ± 1.9) x 10“^. 
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